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[i]  Observations  of  the  electron  density  ne  based  on  measurement  of  the  upper  hybrid 
resonance  frequency  by  the  Polar  spacecraft  Plasma  Wave  Instrument  (PWI)  are  available 
for  March  1996  to  September  1997,  during  which  time  the  Polar  orbit  sampled  all 
MLT  values  three  times.  In  a  previous  study,  we  modeled  the  electron  density  dependence 
along  field  lines  as  ne  =  ne0(Rmax/R)a,  where  ne0  is  the  equatorial  electron  density,  Rmax  « 

LRe  is  the  maximum  geocentric  radius  R  to  any  point  on  the  field  line,  and  ot  =  amodei  = 

8.0  -  3.0  log10«eo  +  0.28(logio«eo)2  -  0.43(Rmax/RE),  for  all  categories  of  plasma 
(plasmasphere  and  plasmatrough).  (In  the  formula  for  otmodei,  ne0  is  expressed  in  cm-3.) 

Here,  we  illustrate  the  field  line  dependence  using  several  example  events.  We  show  that 
the  plasmapause  is  much  more  evident  on  the  large  radius  portion  of  the  orbit  and  that 
at  R  ~  2  RE  the  electron  density  tends  to  level  out  at  large  Rmax  to  a  constant  value 
~100  cm-3.  We  also  present  an  example  of  plasmaspheric  plasma  extending  out  to  at  least 
L  ~  9  on  the  dawnside  during  particularly  calm  geomagnetic  conditions  (as  indicated  by 
low  Kp).  Then  we  present  the  average  equatorial  profiles  of  ne0  versus  Rmax  for 
plasmasphere  and  plasmatrough.  Our  average  plasmasphere  profile  is  found  to  have  values 
intermediate  between  those  based  on  the  models  of  Carpenter  and  Anderson  and  Sheeley 
et  al.  The  plasmatrough  equatorial  density  ne0  scales  with  respect  to  Rmax  like  R^l*,  but  in 
the  region  for  which  our  plasmatrough  data  is  most  reliable  (L  <  6),  it  is  well  fit  by 
the  Rmax°  scaling  of  Sheeley  et  al.  or  the  Rmax5  scaling  of  Carpenter  and  Anderson.  We 
present  a  simple  interpretation  for  the  field  line  dependence  of  the  density.  For  large  neo» 
such  as  occurs  in  the  plasmasphere,  ct  is  close  to  zero  on  average  (implying  that  ne  is 
roughly  constant  along  field  lines).  When  ne0  decreases,  so  does  ne  at  R  =  2  RE,  but  the 
value  there  does  not  decrease  much  below  100  cm-3.  (It  is  unclear  if  this  value  is  an 
absolute  lower  density  limit  because  most  often  the  upper  hybrid  resonance  emission 
disappears  at  R  ~  2  RE  because  fjjfce  <  1,  where  fp  cc^fnt  is  the  plasma  frequency  and  fce  is 
the  electron  cyclotron  frequency.)  Finally,  we  examined  the  dependence  of  ot  and  the 
density  at  the  equator  and  at  R  ~  2  RE  on  the  average  (Kp)  (Kp  averaged  with  a  3-day 
timescale).  There  is  no  clear  dependence  of  the  average  ot  —  amodei  on  (Kp)  or  on  MLT.  In 
the  plasmasphere,  ne0  decreases  with  respect  to  increasing  (Kp).  Index  Terms:  2768 
Magnetospheric  Physics:  Plasmasphere;  2730  Magnetospheric  Physics:  Magnetosphere — inner;  2740 
Magnetospheric  Physics:  Magnetospheric  configuration  and  dynamics;  2794  Magnetospheric  Physics: 

Instalments  and  techniques;  KEYWORDS:  magnetospheric  electron  density,  empirical  density  model, 
plasmasphere 

Citation:  Denton,  R.  E.,  J.  D.  Menietti,  J.  Goldstein,  S.  L.  Young,  and  R.  R.  Anderson  (2004),  Electron  density  in  the  magnetosphere, 
J.  Geophys.  Res.,  109,  A09215,  doi:10.1029/2003JA010245. 


1.  Introduction 

[2]  Modeling  the  electron  density  ne  is  important  because 
ne  can  serve  as  a  proxy  for  the  mass  density  p,  which  affects 
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the  timescale  of  the  magnetospheric  response  to  solar  wind 
perturbations  and  determines  the  properties  of  ULF  waves. 
More  directly,  the  electron  density  can  affect  whistler  waves 
that  scatter  radiation  belt  electrons  and  can  affect  radio 
communications  at  ionospheric  altitudes. 

[3]  The  seminal  work  on  the  distribution  of  the  equatorial 
electron  density  «e0  is  that  of  Carpenter  and  Anderson 
[1992].  More  recently,  models  for  ne 0  have  been  developed 
by  Gallagher  et  al  [2000]  and  Sheeley  et  al  [2001].  The 
latitudinal  density  dependence  along  field  lines  is  less  well 
known.  Methods  used  to  infer  the  latitudinal  dependence  of 
ne  along  magnetic  field  lines  include  in  situ  spacecraft 
observations  and  passive  remote  sensing  with  whistler 
waves.  Another  recent  technique  for  determining  the  field 
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line  distribution  of  ne  is  active  remote  sensing  using  radio 
waves  [Reinisch  et  al ,  2001,  2004].  This  technique  has  the 
advantage  that  instantaneous  profiles  of  «e  along  the  mag¬ 
netic  field  line  can  be  obtained,  but  so  far  the  technique  has 
been  limited  to  cases  where  the  equatorial  electron  density 
nc o  is  large  (several  hundred  cm~3).  Ultra  low  frequency 
(ULF)  toroidal  AlfVen  frequencies  have  been  used  to  infer 
the  field  line  dependence  of  p  (see  references  by  Menk  et  al 
[1999]  and  Denton  et  al.  [2001]  and  also  work  by  Denton  et 
al.  [2004]  and  Takahashi  et  al.  [2004]). 

[4]  Three  previous  studies  have  used  the  Polar  spacecraft 
plasma  wave  data  to  infer  the  typical  dependence  of  electron 
density  ne  along  field  lines.  Denton  et  al  [2002a]  used  a 
method  which  did  not  initially  assume  any  particular 
functional  form  for  the  parallel  dependence  but  then  found 
they  could  fit  their  results  to  the  power  law  form 

nc  =  wco(^max/^)  >  (1) 

where  Rmax  is  the  maximum  geocentric  radius  R  to  any 
point  on  the  field  line  (=  LRE  for  a  dipole  field),  and  «e0  is 
the  value  of  nc  at  R  =  Rmax.  (For  our  data  set,  R  «  i?max 
occurs  at  the  position  of  minimum  magnetic  field  B0 ,  which 
we  call  the  magnetic  equator.)  For  a  dipole  magnetic  field, 
this  form  becomes  «e  =  ne 0  (cos(X))"2a,  which  is  quite 
similar  to  the  form  recently  used  by  Reinisch  et  al  [2004] 
and  Huang  et  al.  [2004],  ne  =  nc0  (cos[(tt/2)  (aX/Xinv)])“p, 
where  Xinv  is  the  invariant  (Earth’s  surface)  value  of  the 
latitude  X.  While  we  do  not  wish  to  claim  that  equation  (1) 
is  the  optimal  functional  form  for  describing  parallel  density 
dependence,  for  a  one  parameter  fit  (a)  it  appears  to  do 
quite  well  [Denton  et  al,  2002a],  especially  considering  the 
large  spread  in  the  parallel  dependence  of  the  data.  Denton 
et  al.  [2002a]  listed  a  at  a  number  of  Rmax  values,  finding, 
for  instance,  a  =  0.8  ±  1.2  at  Rmax  =  4.4  RE  and  a  =  2.1  ± 
1.4  at  Rmax  =  7  Re.  Goldstein  et  al  [2001],  using  a  method 
similar  to  that  of  the  present  paper,  found  an  average  value 
of  a  =  0.37  ±  0.8  for  their  plasmasphere  data  ( ne  > 
100  cm-3)  and  an  average  value  of  a  =  1.7  ±  1.1  for  the 
plasmatrough  («e  <  100  cm"”3;  they  discarded  outlying 
points  before  computing  the  errors).  Denton  et  al.  [2002b] 
used  a  much  larger  data  set  than  that  of  Goldstein  et  al 
[2001].  They  mapped  the  magnetic  field  from  the  spacecraft 
position  at  the  time  of  each  measurement  to  determine  Rmax 
more  accurately.  Then  they  modeled  the  statistical  average 
of  ot  as  a  function  of  the  equatorial  density  ne0  and  Rmax  to 
find  a  single  formula  for  a  for  all  types  of  plasma 
(plasmasphere  and  plasmatrough) 

O^model  =  +  O 

<*»„  =6.0-3.0Ioglon<!o  +  0.28(log10«co)2,  (2) 

°^»oi  =  2-0  —  0A3(Rmm/R£), 

for  2.5  RE  <  Rmax  —  8.5  RE,  2  RE  <  R  5;  Rmaxj  2nd  2  <  ncQ 
<  1500  cm-3.  The  average  error  between  the  observed 

a  values  and  am0dci,  ^/((«  -  cw*)’)  =  0.65.  Denton  et  al. 
[2002b]  stated  that  there  was  no  remaining  dependence  of 
the  average  a  on  MLT  or  Kp  but  did  not  demonstrate  this 
result  due  to  space  limitations. 

[s]  Here,  after  reviewing  the  basic  method  of  analysis 
(section  2),  we  discuss  some  example  events  (section  3).  We 


then  present  a  number  of  average  quantities,  such  as  the 
average  /?max-dependent  profiles  of  n^,  ne  at  R  ~  2  RB  and 
a  (section  4).  In  section  5  we  examine  the  dependence  of  a, 
nc0,  and  nc  at  R  ~  2  RE  on  the  average  (Kp)  (Kp  averaged 
with  a  3-day  timescale).  We  summarize  our  findings  in 
section  6. 


2.  Polar  Electron  Density  Data  and  Method  of 
Analysis 

[6]  The  electron  density  values  used  in  this  paper  are 
obtained  using  the  Polar  Plasma  Wave  Instrument  (PWI) 
[Gurnett  et  al.,  1995],  The  electron  number  density  can 
be  determined  from  noise  emission  which  has  an  upper 
edge  in  frequency  at  the  upper  hybrid  resonance  (UHR) 
frequency  [Goldstein  et  al,  2001],  For  each  data  point,  a 
field  line  mapping  program  was  used  to  map  the  space¬ 
craft  location  to  the  position  along  the  field  line  with 
maximum  radius,  Rmax.  A  Tsyganenko  magnetic  field 
model  was  used  [Tsyganenko,  1995]  as  described  by 
Denton  et  al.  [2002a],  In  cases  for  which  solar  wind 
data  were  available  (92%  of  the  data),  we  used  the 
Tsyganenko  1996  [Tsyganenko,  1995]  magnetic  field 
model  coupled  with  the  International  Geomagnetic  Refer¬ 
ence  Field  (IGRF)  inner  magnetic  field  model  [I AG  A 
Division  V  Working  Group  8,  1991].  Where  the  solar  wind 
data  were  not  available,  we  substituted  the  Tsyganenko 
1989c  model  [Tsyganenko,  1989]  for  the  outer  field. 

[7]  Figure  1  shows  the  path  of  the  Polar  spacecraft  on 
3  February  1997,  projected  onto  a  meridional  plane  (the  Z 
coordinate  is  the  same  as  the  SM  coordinate,  and  X  is  the 
other  meridional  coordinate).  Owing  to  the  nature  of  the 
Polar  orbit,  the  trajectory  of  the  spacecraft  crosses  a 
particular  Rmm  value  (i.e.,  L  shell)  at  two  different  points, 
such  as  “1”  and  “2”  in  the  figure.  At  these  points  we 
can  calculate  (from  equation  (1)) 


Iog(wc2/nci) 
log (/?,//t2)  ’ 


(3) 


where  the  subscripts  1  and  2  indicate  data  on  the  inner  and 
outer  segments  of  the  Polar  trajectory.  (Errors  in  this 
procedure  are  discussed  by  Denton  et  al.  [2002b].) 


3.  Example  Events 

3.1.  The  3  February  1997  Event  With  a  Clear 
Plasmapause 

[s]  Figure  2a  shows  the  electron  density  ne  as  determined 
by  the  PWI  on  3  February  1997, 1322-1453  UT  [Denton  et 
al.,  2002b],  In  Figure  2a,  nt  is  plotted  versus  Rmax  for  the 
outer  (thin  curve)  and  inner  (bold  curve)  portions  of  the 
orbit.  (In  the  example  plots,  the  thicker  curve  corresponds  to 
the  high  density  portion  of  the  orbit.)  Values  of  the  radius  R 
for  both  portions  are  plotted  in  Figure  2b.  The  resolution  of 
the  PWI  data  points  is  ~0. 1  RE,  as  described  by  Denton  et  al. 
[2002b],  While  it  is  clear  that  there  is  some  azimuthal  or 
temporal  dependence  (the  plasmapause  position  is  not  ex¬ 
actly  the  same  for  the  large  and  small  radius  portions  of  the 
orbit),  in  this  particular  case  that  dependence  appears  to  be 
small.  Furthermore,  the  smooth  variation  of  density  on  both 
portions  of  the  orbit  gives  evidence  that  in  this  case  there  is 
not  a  great  amount  of  structure  in  the  azimuthal  direction. 
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Figure  1.  The  path  of  the  Polar  spacecraft  on  3  February 
1997  projected  onto  a  meridional  plane.  From  the  values  of 
radius  R  and  electron  density  we  at  points  1  and  2, 
equation  (1)  can  be  solved  for  a  and  nc0- 

[9]  The  a  values  plotted  in  Figure  2c  are  determined  at  a 
number  of  Rnmx  values  spaced  approximately  0.1  RE  apart 
and  within  the  range  of  Rmax  represented  by  the  data  on  both 
inner  and  outer  radius  segments  (here  Rmax  =  3.0-6. 0  RE). 
Note  that  the  value  of  a  is  small  within  the  plasmasphere 
where  the  density  on  the  inner  and  outer  radius  portions  of 
the  orbit  are  nearly  the  same.  A  value  of  a  —  0  corresponds 
to  a  flat,  or  constant,  density  dependence  (equation  (1)).  For 
this  particular  event,  a  is  also  small  just  outside  the 
plasmapause,  but  it  increases  farther  out  in  the  plasmatrough 
where  the  density  on  the  inner  and  outer  radius  portions  of 
the  orbit  differ. 

[10]  The  variation  of  density  along  the  field  line  can  also 
be  described  in  terms  of  a  scale  length.  A  physically 
meaningful  definition  of  a  scale  length  might  be  L ^  = 
y/2ntQ /{dP-nJdP )  |/==0,  where  /  is  the  length  along  the  field 
line  defined  to  be  zero  at  the  magnetic  equator.  Assuming 
dnJdl\i=o  =  0  at  the  magnetic  equator  (implied  by  equation 
(1)),  L ^  is  the  distance  along  the  field  line  from  the  equator 
at  which  the  density  is  twice  the  equatorial  value.  In  a  dipole 
field,  -  (LRE)ly/ a.  Based  on  this  fact,  we  define  (for 
arbitrary  magnetic  field) 

R<x  —  Rmax/  y/o.  (4) 


and  plot  values  of  £<*  in  Figure  2d.  As  was  the  case  for  a, 
values  of  La  are  not  meaningful  in  the  vicinity  of  the 
plasmapause  where  the  gradients  in  ne  with  respect  to  Rmax 
are  large.  Outside  the  plasmapause  ( Rmax  £  3.6  RE),  ZQ  is 
roughly  constant  RE  in  this  case. 

[11]  One  notable  feature  of  Figure  2a  is  that  the  density 
on  the  low  radius  portion  of  the  orbit  (thick  curve 
corresponding  to  R  ~  2  RE )  levels  off  at  ~100  cm-3  at 
large  Rmax.  This  is  a  common  feature  that  we  see  in  many 
Polar  orbits  and  indicates  that  at  large  Rmsx  the  density  is  a 
stronger  function  of  radius  than  of  L  shell.  Values  of  ne  ~ 
100  cm-3  at  R  ~  2  RE  with  predominantly  radial  depen¬ 
dence  are  also  found  in  the  polar  cap  [Persoon  et  al,  1983; 
Gallagher  et  al ,  2000;  Nsumei  et  al ,  2003].  While  at  this 
altitude  (R  ~  2  RE),  the  density  at  large  Rmax  (but  still  on 
closed  field  lines)  becomes  more  polar-cap-like,  the  density 
at  the  same  Rmax  values  but  in  the  vicinity  of  the  magnetic 
equator  has  different  properties  than  the  polar  cap  density. 
The  differences  are  the  Rmax  dependence  (thin  curve  in 
Figure  2a),  the  fact  that  the  radial  dependence  along  the 
field  lines  is  not  so  steep  [Denton  et  al ,  2002b;  Nsumei  et 
al,  2003],  and  the  fact  that  ne0  decreases  or  is  at  least  flat 
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Figure  2.  As  a  function  of  Rmax,  (a)  the  electron  density  «e 
as  determined  from  the  PWI,  (b)  the  Polar  geocentric 
radius  R,  (c)  the  power  law  index  a  (equation  (1)),  and 
(d)  the  parallel  scale  length  (equation  (4))  as  determined 
from  Polar  data  measured  on  3  February  1997,  1322— 
1453  UT.  In  Figures  2a  and  2b,  the  thin  (bold)  curve 
represents  the  profile  of  ne  determined  at  the  large  (small) 
radius  portion  of  the  orbit. 
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Figure  3.  Same  as  Figure  2,  except  for  data  collected  on 
13  January  1997,  2140-2320  UT. 

with  respect  to  Kp  (sections  5.2  and  5.3)  rather  than 
increasing  with  respect  to  Kp  [Nsumei  et  al.,  2003].  All  of 
these  differences  relate  to  trapping  of  particles  in  the 
plasmatrough  and  consequent  buildup  of  ne0  on  closed  field 
lines.  Even  at  R  ~  2  RE  the  Kp  dependence  appears  to  be 
somewhat  different  as  shown  in  section  5.2.  The  plasma- 
trough  might  then  be  considered  a  transition  region  sharing 
some  properties  of  the  polar  cap  and  the  plasmasphere. 

3.2.  Events  for  Which  the  Plasmapause  Is  Not  So 
Apparent  on  the  Low-Radius  Portion  of  the  Orbit 

[n]  Figure  3  shows  an  orbit  with  a  plasmapause  that  is 
less  steep.  Though  the  inner  and  outer  portions  of  the  orbit 
have  significantly  different  R  (Figure  3b),  the  density 
structure  on  the  two  portions  of  the  orbit  is  quite  similar 
(Figure  3a).  However,  the  density  values  for  the  inner  curve 
(bold  curve)  do  not  drop  to  the  low  values  of  the  outer  curve 
(thin  curve  with  ne  ~  4  cm"3  at  Rmm  >  5  RE).  In  this  case, 
the  density  on  the  inner  curve  decreases  to  ~50  cm-3  at 
•Kmax  =  5.2  Re.  Correspondingly,  the  value  of  a  is  close  to 
zero  at  low  J?max5i  4  Rp  but  is  larger  outside  the  plasma¬ 
pause  at  Rm3yk,  5  RE- 

[13]  Figure  4  shows  an  orbit  where  the  density  just  inside 
the  plasmapause  is  not  very  far  above  100  cm"3.  In  this 
case,  the  plasmapause  is  evident  on  the  high-radius  portion 
of  the  trajectory  (thin  curve)  but  not  on  the  low-radius 
portion  of  the  trajectoiy  (thick  curve).  Again,  the  density 
does  not  decrease  very  much  below  about  100  cm"3  for  the 


low-radius  portion  of  the  orbit  with  R~2RE  (at  least  in  the 
vicinity  of  the  plasmapause  at  /?max  ~  5.5  RE;  the  behavior 
at  the  largest  f?max  values  >7  Rp  is  not  clear;  see  section  3.4). 
While  ne  often  levels  off  to  ~100  cm"3  on  the  inner  portion 
of  the  Polar  orbit  (with  R  ~  2  RE)  at  large  Rmm,  ne  can 
become  quite  large  on  the  inner  portion  of  the  orbit  at  low 
^max  where  the  flux  tubes  are  filled,  that  is,  in  the  plasma- 
sphere  (Figures  2a,  3a,  and  4a). 

3.3.  The  19  December  1996  Event  With  No 
Plasmapause 

[m]  Figure  5  shows  data  collected  on  19  December 
1996,  2154-2440  UT.  This  event  falls  in  the  “gradually 
decreasing”  category  of  Denton  et  al.  [2002b]  (see  Ap¬ 
pendix  A).  The  density  gradually  decreases  on  both  inner 
and  outer  portions  of  the  orbit  so  that  no  clear  plasma- 
pause  is  evident.  Another  possible  interpretation  of  this 
data  (as  pointed  out  by  a  referee)  is  that  the  plasmapause  is 
very  wide,  ending  at  ~/?max  ~  7  Rp.  However,  our  opinion 
is  that  there  is  not  much  difference  between  a  very  wide 
plasmapause  and  no  plasmapause  at  all.  It  is  not  clear  how 
to  distinguish  a  very  wide  plasmapause  from  the  decrease 
in  density  with  respect  to  Rmm  within  the  plasmasphere  or 
plasmatrough.  This  event  includes  measurements  out  to 
■ftmax  =  9.8  Re  based  on  the  Tsyganenko  1996  magnetic 
field  model  [Tsyganenko,  1995]  (out  to  L  =  8.7  based  on 
the  dipole  field  model).  The  reason  the  spacecraft  observes 


-^max  (*E) 


Figure  4.  Same  as  Figure  2,  except  for  data  collected  on 
3  October  1996,  0401-0645  UT. 
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Figure  5.  Same  as  Figure  2,  except  for  data  collected  on 
19  December  1996,  2154-2440  UT. 

such  large  Rmax  values  on  the  outer  portion  of  the  orbit  is 
that  in  this  case  the  magnetic  latitude  MLAT  there  is  equal 
to  33.4°. 

[is]  It  appears  from  Figure  5  that  the  plasmasphere,  at 
least  in  the  sense  of  the  “gradually  decreasing”  plasma 
category  (Appendix  A),  can  extend  out  to  large  Rmax  (Rmax 
~  9-10  Re).  Some  observations  of  large  density  at  large 
Rmax  correspond  to  plasma  plumes  [Sandel  et  al ,  2001], 
which  usually  occur  on  the  duskside.  However,  the  obser¬ 
vation  shown  in  Figure  5  occurred  at  magnetic  local  time 
MLT  =  7  hours.  Furthermore,  similar  profiles  of  density 
with  respect  to  Rmax  were  observed  at  the  same  MLT  on 
the  next  two  Polar  passes  and  at  dusk  MLT  ~20  hours, 
before  and  after  the  time  corresponding  to  the  event,  as 
described  in  Appendix  B.  This  extended  plasmasphere 
occurred  during  an  extended  period  of  low  geomagnetic 
activity  as  indicated  by  low  Kp  (see  Appendix  B),  a 
condition  which  normally  allows  the  plasmasphere  to 
expand.  It  is  possible  that  the  large  densities  are  confined 
to  azimuthally  localized  structures  [Sandel  et  al ,  2001]. 
On  the  other  hand,  the  number  of  observations  at  different 
times  and  at  different  MLT  values  (Appendix  B)  may 
indicate  that  the  plasma  density  has  increased  to  plasma- 
spheric  levels  throughout  the  magnetosphere  Rmax  < 
10  RE.  At  the  least,  the  high-density  region  (or  regions) 
must  fill  a  significant  fraction  of  the  magnetosphere.  In 
any  case,  regardless  of  whether  the  observed  density  is 
azimuthally  localized  or  not  this  example  clearly  shows 


that  there  are  cases  for  which  there  is  no  clear  plasma- 
pause,  at  least  along  a  radial  cut,  and  in  this  sense  the 
plasma  at  these  large  Rmax  values  is  plasmasphere-like. 

3.4.  Loss  of  Data  at  Large  Rmax 

[  16]  One  common  feature  we  see  in  all  the  example 
events  (Figures  2-5)  is  the  cutoff  of  the  data  at  large  Rm ax 
on  the  low-radius  portion  of  the  orbit;  in  each  case,  the  data 
on  the  large  radius  portion  of  the  orbit  (thin  curves  in 
Figures  2a,  3a,  4a,  and  5a)  continues  to  larger  Rmax  values 
than  the  data  on  the  low  radius  portion  (thick  curves  in 
Figures  2a,  3a,  4a,  and  5a).  For  instance,  in  Figure  2a,  the 
data  on  the  low-radius  portion  of  the  orbit  (thick  curve)  ends 
at  Rmax  =  6.0  RE,  whereas  the  data  on  the  high-radius 
portion  (thin  curve)  continues  to  beyond  the  highest  Rmax 
value  plotted  ( Rmax  =  6.4  RE).  We  have  investigated  this 
data  cutoff  for  the  four  events.  In  each  case,  it  occurs  when 
fplfce  ~  L  where  fp  is  the  plasma  frequency  oc^/n^  and  fce  is 
the  electron  cyclotron  frequency  proportional  to  the  mag¬ 
netic  field  B.  The  existence  of  the  upper  hybrid  noise  band 
generally  requires  fplfce  1  [Tataronis  and  Crawford,  1970]. 

[17]  In  the  case  of  Figure  2  at  MLT  ~15  hours,  the  data 
cutoff  at  low  radius,  R  =  2.0  RB,  occurs  at  Rmax  =  6  RB.  The 
density  could  be  lower  outside  this  location  (leading  to 
lower  fpffceOty/fQ,  or  the  reduction  in  fp/fce  could  be  due  to 
the  increasing  fceccB .  The  particle  data  (Polar  TIMAS)  and 
wave  properties  (PWI)  indicate  that  Polar  continues  to  be  in 
the  plasmatrough  for  some  distance  outside  of  Rmsx  =  6.0  RE. 
We  simply  do  not  know  what  the  density  is  during  this  time. 
AtL  =  9.6  (based  on  a  dipole  field),  Polar  enters  the  auroral 
region,  as  indicated  by  the  sudden  appearance  of  enhanced 
broadband  electrostatic  waves  and  electromagnetic  auroral 
hiss  associated  with  this  regions  [Gumett  and  Inan ,  1988; 
Lin  et  al,  1984]. 

[is]  In  the  case  of  Figure  3,  at  MLT  ~5.6  hours,  the 
disappearance  of  the  upper  hybrid  emission  at  low  radius, 
R  =  2.4  Re,  occurs  at  Rmax  =  5.2  i?E.  Again,  we  do  not  know 
what  the  density  is  at  larger  Rmax.  In  this  case,  however,  the 
disappearance  of  the  upper  hybrid  emission  coincides  with 
entry  into  the  auroral  region  or  cusp,  as  noted  by  a  sudden 
change  in  the  plasma  wave  spectrum  as  described  in  the  last 
paragraph. 

[19]  The  data  associated  with  Figure  4  is  similar  to  that 
for  Figure  2.  The  largest  Rmax  value  for  which  we  have  a 
density  measurement  is  at  Rmax  =  7.2  RE  (R  =  2.3  RE  and 
MLT  ~  11  hours).  There  is  a  transition  to  auroral  zone 
properties  at  L  ~  11  (based  on  a  dipole  field).  There  could 
be  a  region  of  the  plasmatrough  in  between  the  two  regions 
for  which  we  cannot  measure  the  density.  However,  unlike 
the  case  of  Figure  2a,  in  Figure  4a  the  density  is  apparently 
decreasing  at  the  time  that  fp/fce  ~  1  and  the  upper  hybrid 
emission  disappears.  Therefore  the  plasmatrough  density 
may  very  well  decrease  below  100  cm-3. 

[20]  In  the  case  of  Figure  5,  corresponding  to  MLT  =  7.5 
hours,  the  low-radius  data  cutoff  is  at  Rmax  =  9.6  RE  (at  R  ~ 
2.2).  This  case  is  like  that  of  Figure  3.  The  disappearance  of 
the  upper  hybrid  emission  coincides  with  entry  into  the 
auroral  region  as  noted  by  a  sudden  change  in  the  plasma 
wave  spectrum. 

[21]  In  summary,  the  low-radius  data  cutoff  was  appar¬ 
ently  associated  with  a  drop  in  fpffce  below  unity.  For  two 
out  of  four  of  the  events,  our  low-radius  data  cutoff  occurred 
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Figure  6.  Equatorial  electron  density  «e0  versus  Rm ax.  For 
the  plasmasphere  (thick  curves),  we  plot  the  log  average  «e0 
based  on  the  Polar  data  (thick  solid  curve),  the  saturated 
plasmasphere  model  of  Carpenter  and  Anderson  [1992] 
(thick  dotted  curve),  the  plasmasphere  model  of  Sheeley  et 
al.  [2001]  (thick  short  dashed  curve),  and  the  plasmasphere 
model  of  Gallagher  et  al  [2000]  (thick  long  dashed  curve). 
For  the  plasmatrough  (thin  curves),  we  plot  the  log  average 
«e0  based  on  the  Polar  data  (thin  solid  curve),  the  best  power 
law  fit  to  our  data  (thin  dotted  curve),  and  the  plasmatrough 
model  of  Sheeley  et  al  [2001]  (thin  dashed  curve).  The  dot- 
dashed  curves  are  the  log  average  ne0  based  on  the  Polar 
data  but  using  the  Sheeley  et  al  [2001]  density  criterion 
(equation  (14))  to  distinguish  the  data  for  the  plasmasphere 
(thick  dot-dashed  curve)  from  the  plasmatrough  (thin  dot- 
dashed  curve). 


in  association  with  passage  into  the  auroral  zone.  For  the 
other  two  events,  there  may  have  been  a  region  of  the 
plasmatrough  with  decreased  density  outside  of  the  region 
for  which  we  were  able  to  measure  the  density  based  on 
the  presence  of  the  upper  hybrid  resonance  band.  In  one 
of  these  cases  (Figure  4),  the  density  appears  to  be 
dropping  when  the  upper  hybrid  emission  ceases.  In  the 
other  case  (Figure  2),  there  could  be  a  decrease  in  density, 
or  the  change  in  f/fce  could  be  explained  by  increasing  B 
at  low  R . 

4.  Rmax  Dependence 
4.1.  Rmax  Dependence  of  we0 

[22]  Figure  6  shows  the  Rmax  dependence  of  the  log 
average  equatorial  electron  density  nc0  for  plasmasphere 
(thick  solid  curve)  and  the  plasmatrough  (thin  solid  curve). 
The  categories  of  plasma  were  defined  by  visual  inspection 
as  described  in  Appendix  A.  For  the  purposes  of  this  paper, 
the  plasmasphere  data  set  includes  the  “gradually  decreas¬ 
ing”  data  set  of  Denton  et  al  [2002b]  (see  section  3.3  and 
Appendix  A);  the  “gradually  decreasing”  plasma  is  as¬ 
sumed  to  be  extended  plasmasphere. 


[23]  The  thick  dotted  curve  in  Figure  6  is  Carpenter  and 
Anderson's  [1992]  saturated  plasmasphere  density 

logjo  ("co/cirT3)  =  fl]  +a2^max 

+  (0.00127*sun  -  0.0635)e-^-2)/L5,  (5) 

where  RSuu  is  the  average  sunspot  number  ~13  over  the 
period  of  our  data,  we  have  averaged  over  their  sinoida! 
local  time  dependence,  Carpenter  and  Anderson’s  values  for 
a  1  and  a2  are 


=  3.90} 

=  —0.315, 

and  we  have  substituted 

Rmax=Rmax/RE  (7) 

(~  L  for  a  dipole  field)  for  their  L .  (The  Carpenter  and 
Anderson  study  used  spacecraft  data  with  MLAT  <  30°,  and 
the  distinction  between  L  based  on  a  dipole  field  versus  our 
^max  is  likely  to  be  less  important  than  for  our  data  set  from 
the  Polar  spacecraft.)  The  short  dashed  curve  in  Figure  6  is 
the  plasmasphere  model  of  Sheeley  et  al  [2001] 

*e0  =  ( 1 390cm“3)  (3/ ^max  )4  83 , 

±(440cm-3)(3/^max)36°  (8) 


for  3  <  *max  <  7. 

[24]  Our  equatorial  average  plasmasphere  ne0  values 
(thick  solid  curve)  lie  between  the  nc0  values  of  Carpenter 
and  Anderson  (thick  dotted  curve  based  on  equation  (5) 
with  equation  (6))  and  those  of  Sheeley  et  al.  (thick  short- 
dashed  curve  based  on  equation  (8)).  Our  plasmasphere  data 
can  be  best  fit  using  equation  (5)  with 


a\  =3.78  ±0.18, 
a2  =  -0.324  ±  0.038, 


(9) 


for  2.5  <  Rmax  <  8.0.  Our  coefficients  (equation  (9))  agree 
with  Carpenter  and  Anderson’s  values  within  the  uncertain¬ 
ties.  (The  fact  that  the  Polar  ne0  values  are  lower  than  values 
from  the  other  models  at  low  flmax  <  3.5  RE  may  be  related 
to  the  boom  oscillation  problem  on  Polar  [Denton  et  al, 
2002b].)  On  average,  the  Polar  ne0  measurements  are  a 
factor  of  1.9  off  from  our  fit,  equation  (5)  with  equation  (9) 
(varying  in  Rmax  bins  from  a  factor  of  1.4  off  at  flmax  =  2.7 
to  a  factor  of  2.0-2.5  off  at  Rmax  >  4.9). 

[25]  The  long  dashed  curve  in  Figure  6  is  the  plasma¬ 
sphere  model  of  Gallagher  et  al  [2000],  which  is  equation 
(5)  with  a  1  and  a2  given  by 


a™  =  5.3, 
a™  =  -0.79. 


(10) 


(Actually,  Gallagher  et  al.’s  equation  for  «e0  has  extra  terms 
other  than  equation  (5),  but  equation  (5)  captures  the  main 
dependence  of  their  model.)  There  is  clearly  a  significant 
difference  between  the  plasmasphere  model  of  Gallagher  et 
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al  [2000]  and  the  other  models.  Gallagher  et  al.’s  model 
(equation  (5)  with  equation  (10))  is  based  on  ion  particle 
data  and  might  not  have  measured  a  significant  fraction  of 
the  electron  density  at  large  i?max.  On  the  other  hand,  the 
data  used  (DE-1  RIMS)  was  calibrated  against  measure¬ 
ments  from  the  DE  1  Plasma  Wave  Instrument  (D. 
Gallagher,  private  communication,  2004).  Another  possible 
explanation  for  the  difference  relates  to  the  method  by 
which  the  data  was  screened.  Gallagher  et  al.’s  data  set  is 
based  on  low  Kp,  with  the  3-day  weighted  average  (section 
5.1)  <1.3.  They  may  have  included  data  in  their  plasma- 
spheric  data  set  that  we  would  have  classified  as 
plasmatrough  based  on  visual  inspection  (see  Appendix  A). 

[26]  As  was  already  mentioned,  the  thin  solid  curve  in 
Figure  6  is  the  log  average  plasmatrough  density  based  on 
the  Polar  spacecraft  data.  Carpenter  and  Anderson's  [1992] 
plasmatrough  model  uses  a  power  law  form 

ne0  =  ai(Rmix)a\  (11) 

with  a4  =  -4.5  and  a3  specified  by  a  condition  matching 
plasmasphere  and  plasmatrough  density  (we  have  substi¬ 
tuted  i?max  for  their  L).  Our  plasmatrough  data  is  best  fit 
with  a4  =  — 3.45, 

log,o  («eo/cm-3)  =  (3.77  ±  0.31)  +  (-3.45  ±  0.41)  log)0 

(12) 

fit  for  2.9  <  ^max  <  7.5,  and  plotted  as  the  thin  dotted  curve 
in  Figure  6.  On  average,  the  Polar  density  measurements  are 
a  factor  of  2.1  off  from  this  model.  Sheeley  et  al  [2001] 
found  a4  =  —4.0.  The  value  a4  =  —4.0  is  what  one  would 
expect  if  the  flux  tube  volume  (per  magnetic  flux)  oc  L(B  ~ 
jR*  ax  and  the  flux  tube  particle  content  (per  magnetic  flux)  is 
constant  across  the  plasmatrough  (Appendix  C).  Their 
model  for  plasmatrough  density  is 

He 0  =  124  cm-3(3/ijma*)4  0, 

±  78  cm~3(3/,ftmax)4'72,  (13) 

for  3  <  Rmax  <  7,  where  we  have  averaged  over  their 
sinoidal  local  time  dependence  and  have  again  substituted 
Rmsx  for  L.  Their  plasmatrough  ne0  is  plotted  as  the  thin 
dashed  curve  in  Figure  6.  Our  curve  (thin  solid  curve)  lies 
very  close  to  that  of  Sheeley  et  al.  (thin  dashed  curve) 
except  at  large  /?max  where  our  curve  levels  out  to  ~He0  ~ 
7  cm-3  (we  will  comment  on  this  in  the  last  paragraph  of 
this  subsection). 

[27]  Sheeley  et  al *s  [2001]  study  used  plasma  wave  data 
as  does  ours  (their  data  is  found  from  observations  by  the 
CRRES  spacecraft).  One  difference  in  Sheeley  et  al.’s  data 
analysis  is  that  they  distinguished  their  plasmasphere  and 
plasmatrough  data  using  a  criterion  on  density 

He 0  =  (10  cm_3)(6.6/^max)4  (14) 

(correcting  a  misprint  in  their  paper).  Observations  with 
density  above  (below)  this  value  were  assumed  to  be 
plasmasphere  (plasmatrough).  We  separated  plasmasphere 
and  plasmatrough  data  for  individual  events  using  a  more 
complicated  procedure  (Appendix  A),  similar  to  that  of 
Carpenter  and  Anderson  [1992].  To  test  whether  the 


difference  in  results  is  due  to  this  difference  in  method,  we 
also  separated  plasmasphere  and  plasmatrough  data  (from  all 
types  of  plasma)  using  the  criterion  in  equation  (14).  The 
density  values  found  using  this  method  are  plotted  as  dot- 
dashed  curves  in  Figure  6.  The  thick  dot-dashed  curve 
(plasmasphere  chosen  with  equation  (14))  is  close  to  and 
actually  higher  than  the  thick  solid  curve  (our  plasmasphere 
data  set).  Thus  this  difference  in  method  cannot  explain  why 
the  Sheeley  et  al.  «e 0  (thick  short-dashed  curve)  is  lower  than 
ours  (thick  solid  curve). 

[28]  The  thin  dot-dashed  curve  in  Figure  6  (plasmatrough 
determined  by  equation  (14))  is  somewhat  lower  than  our 
thin  solid  curve  (log  average  for  plasmatrough  data  set).  The 
use  of  equation  (14)  to  choose  small  densities  results  in  a 
lower  average  density.  This  could  help  explain  the  differ¬ 
ence  in  density  values  between  the  Sheeley  et  al.  model 
(thin  dashed  curve)  and  our  log  average  plasmatrough 
density  (thin  solid  curve)  at  some  values  of  7?max,  such  as 
i?max  >  5  Re .  On  the  other  hand,  the  Polar  PWI  signal  often 
drops  out  when  the  density  is  very  low  at  large  i?max 
(particularly  when  R  is  low;  see  section  3.4),  and  this  may 
be  the  cause  of  the  leveling  out  of  our  log  average 
plasmatrough  density  to  ~weo  ^  7  cm^3  at  2?max  =  7-7.5 
Rb .  If  we  use  our  plasmatrough  data  only  out  to  Rmax  =  6.4 
7?e,  our  best  fit  is 

log io  (W«n-3)  =  (4.13  ±  0.46)  +  (-4.10  ±  0.74)  log,0 

(15) 

and  now  the  coefficient  a4  =  —4.10  ±  0.74  (see  equation 
(11))  is  equal  to  both  the  Sheeley  et  al.  a4  -  —4.0  and 
Carpenter  and  Anderson  a4  —  —4.5  values  within  the 
uncertainty. 

4.2.  Rmax  Dependence  of  Other  Quantities 

[29]  In  this  subsection,  we  examine  the  i?max  dependence 
of  other  quantities  and  explore  the  result  of  combining  our 
field  line  density  dependence  model  (equation  (1)  with 
equation  (2))  with  the  ^max-dependent  plasmasphere  and 
plasmatrough  «e0  profiles  described  in  section  4.1  (solid 
curves  in  Figure  6).  In  Figure  7  a  number  of  quantities  are 
plotted  versus  I?max.  The  thick  curves  correspond  to  the 
plasmasphere,  while  the  thin  curves  correspond  to  the 
plasmatrough.  The  solid  curves  result  from  an  average. 
The  solid  curves  in  Figure  7a  are  the  same  log  average 
72 e0  values  as  were  plotted  in  Figure  6.  The  solid  curves  in 
Figure  7b  are  the  average  a  values.  The  dashed  curves  all 
relate  to  the  log  average  values  of  ne  computed  for  R  =  2  Re 
(more  precisely,  for  1.8  i?E  <  R  <  2.2  RE).  The  Polar 
spacecraft  perigee  was  at  ~R  -  1.87  (at  the  beginning  of 
1 997),  so  R  -  2  Re  is  essentially  the  lowest  radius  sampled. 
The  dashed  curve  in  Figure  7a  shows  these  log  average 
values.  The  dashed  curves  in  Figure  7b  show  the  values  of 
a  found  from  equation  (3)  using  the  log  average  values  of 
ue  from  Figure  7a  at  the  equator  and  at  R  =  2  RE  (solid  and 
dashed  curves  from  Figure  7a).  The  dotted  curves  all  relate 
to  the  o^odei  model  (equation  (2)).  The  dotted  curves  in 
Figure  7b  are  the  values  of  amodel  with  the  log  average  ne0 
values  (solid  curves  in  Figure  7a)  as  an  input.  The  dotted 
curves  in  Figure  7a  are  the  «e  values  at  R  =  2  RE  based  on 
these  same  am0dei  values  (dotted  curves  in  Figure  7b)  and 
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log  average  equatorial  electron  density  ne 0,  the  dashed 
curves  are  the  log  average  values  of  ne  at  R  =  2  RE,  and  the 
dotted  curves  are  the  values  of  ne  at  R  =  2  RE  assuming  the 
field  line  dependence  (equation  (1))  with  a  =  a^dei 
(equation  (2))  and  using  the  log  average  (solid  curve)  as 
an  input  for  ne0.  In  Figure  7b  the  solid  curves  are  the 
average  o;mot|ci  values,  the  dashed  curves  are  the  values  of  a 
found  from  equation  (3)  using  the  log  average  values  of  ne 
from  (a)  at  the  equator  and  at  R  =  2  RE  (solid  and  dashed 
curves  from  (a)),  and  the  dotted  curves  are  the  amode,  values 
(equation  (2))  as  a  function  of  the  log  average  neo  plotted  in 
Figure  7a  (solid  curves  in  Figure  7a).  In  Figure  7c  the 
curves  are  the  values  of  the  parallel  scale  length  La  found 
from  equation  (4)  using  the  corresponding  a  values  from 
Figure  7b. 


the  log  average  ne0  values  (solid  curves  in  Figure  7a).  The 
parallel  scale  length  values  in  Figure  7c  are  calculated  from 
equation  (4)  using  the  corresponding  a  values  plotted  in 
Figure  7b.  In  other  words,  the  thick  solid  curve  in  Figure  7c 
is  the  value  of  La  based  on  the  thick  solid  curve  from 
Figure  7c,  that  is,  the  La  value  based  on  the  average  a  for 
the  plasmasphere. 

[30]  Figure  7b  shows  that  the  a  values  are  relatively 
small,  roughly  0-1,  in  the  plasmasphere,  while  values 
around  2.5  are  typical  in  the  plasmatrough.  Correspondingly, 


La  is  larger  in  the  plasmasphere,  ~6  RE  for  4  <  rtmax  <  7, 
as  compared  to  the  plasmatrough  with  La  ~  2.5-5  RE. 
( Goldstein  et  al.  [2001]  found  La  ~  5.5  RE  averaging  over 
all  types  of  plasma). 

[31]  As  mentioned  previously,  the  dotted  curves  in 
Figure  7a  are  the  ne  values  at  R  =  2  RE  using  the  amode, 
model  with  the  log  average  n <*,  values  as  an  input.  Both  the 
curves  for  the  plasmasphere  (thick  dotted  curve)  and  plasma- 
trough  (thin  dotted  curve)  approach  values  close  to  100  cm-3 
at  large  Rmax.  Surprisingly,  the  plasmatrough  curve  levels  off 
at  somewhat  higher  density  (~200  cm-3)  than  the  plasma¬ 
sphere  curve  (~100  cm-3).  The  log  average  values  of  ne  at 
•Rmax  =  2  Re  also  approach  a  value  close  to  100  cm-3  at  large 
^max  (dashed  curves  in  Figure  7a).  However,  in  this  case,  the 
curve  for  the  plasmasphere  (thick  dashed  curve  in  Figure  7a) 
is  everywhere  higher  than  the  curve  for  the  plasmatrough 
(thin  dashed  curve  in  Figure  7a),  as  expected. 

[32]  The  difference  in  the  dotted  and  dashed  curves  in 
Figure  7a  demonstrates  the  limitations  of  our  model  as 
expressed  by  equation  (1)  with  equation  (2).  Seeing  as  the 
dotted  curve  in  Figure  7a  is  the  equatorial  density  extrap¬ 
olated  to  R  —  2  Re  and  the  dashed  curve  is  the  actual  density 

R  =  2  Re,  the  fact  that  the  dotted  and  dashed  curves  are 
almost  the  same  for  the  plasmasphere  data  (thick  curves) 
indicates  that  otmoejej  (equation  (2))  can  be  used  to  accurately 
extrapolate  plasmaspheric  equatorial  densities  down  to 
R  -  2  RE.  On  the  other  hand,  for  the  plasmatrough  (thin 
curves  in  Figure  7a),  the  difference  between  the  thin  dotted 
and  dashed  curves  indicates  a  possible  error  on  the  order  of  a 
factor  of  2  at  the  larger  Rmm  values  (5  RE)  when  extrapo¬ 
lating  equatorial  densities  to  R  =  2  RE.  (It  is  not  surprising 
that  the  agreement  for  the  plasmasphere  data  is  better 
considering  that  the  plasmasphere  a  values  are  low  so  that 
the  field  line  dependence  is  weak.)  Correspondingly,  the 
a  values  found  three  different  ways  and  plotted  in  Figure  7b 
agree  fairly  well  for  the  plasmasphere  data,  but  for  the 
plasmatrough,  there  is  disagreement  on  the  order  of  0.5. 

[33]  Based  on  the  observations  described  in  section  3.4, 
we  need  to  qualify  our  results.  It  is  likely  that  sometimes  n[ 
is  lower  than  100  cm-3  at  low  radius  and  large  Rmm  and  we 
are  not  able  to  measure  the  density  because  fp!fce  <  1  so  that 
the  upper  hybrid  noise  band  is  not  observed.  This  does  lead 
to  some  doubt  about  our  conclusion  that  the  density  at  R  ~  2 
Re  does  not  go  significantly  below  100  cm-3.  On  the  other 
hand,  we  often  see  the  density  level  off  to  a  value  near 
100  cm-3  when  the  equatorial  density  is  still  decreasing 
with  respect  to  Rumx  (Figures  2a  and  5a).  Furthermore,  if  we 
were  to  add  observations  with  significantly  lower  density  at 
R  ~  2  Re,  that  would  have  the  effect  of  lowering  the  a 
values  calculated  using  the  average  density  at  the  equator 
and  at  R  =  2  RE  (dashed  curve  in  Figure  7b).  However,  for 
the  plasmatrough  data  (for  which  low  density  at  R  ~  2  RE  is 
more  likely),  these  values  are  already  lower  than  the  average 
a  values  based  on  all  the  measurements,  including  those  at 
larger  R  (thin  solid  curve  in  Figure  7b). 

4.3.  Meridional  Dependence 

[34]  Features  of  the  field  line  dependence  are  shown  in  a 
different  way  in  Figure  8.  Here,  we  plot  in  a  meridional  plane 
contour  lines  for  our  plasmatrough  density  (equation  (12) 
with  equation  (2))  for  the  field  line  dependence.  The  thin 
solid  curves  are  the  density  contours,  the  thick  solid  curves 


8  of  14 


A09215 


DENTON  ET  AL.:  ELECTRON  DENSITY  IN  THE  MAGNETOSPHERE 


A09215 


t 


X(Re) 

Figure  8.  Meridional  plane  showing  contour  lines  (thin 
solid  curves)  of  the  density  given  by  equation  (12)  with 
equation  (2)  for  the  field  line  dependence.  The  thick  solid 
curves  are  dipole  field  lines,  and  the  thick  dashed  curves  are 
circular. 


are  dipole  field  lines,  and  the  thick  dashed  curves  are  circular 
(the  curve  at  a  radius  of  1 RE  is  the  Earth’s  surface).  Note  that 
in  the  region  close  to  the  magnetic  equator  the  density 
contour  that  crosses  the  magnetic  equator  at  low  L  =  3  has 
a  similar  shape  to  the  dipole  field  because  a  is  small  for  the 
large  density  at  this  L  (implying  a  flat  density  along  magnetic 
field  lines).  On  the  other  hand,  the  contour  line  that  crosses 
the  magnetic  equator  at  larger  L  =  7  is  more  circular  in  shape 
(implying  a  density  dependence  that  is  closer  to  a  function  of 
radius).  Note  also  that  the  contour  line  for  L  =  3  has  nearly 
constant  radius  atR~2RE  (position  of  arrow).  This  feature 
corresponds  to  the  leveling  out  of  the  density  with  respect  to 
i?max  at  R  ~  2  RE  (Figure  7a). 

4.4.  Simple  Model  for  a 

[35]  Equation  (2)  describes  the  field  line  variation  of  nc 
(a)  in  terms  of  ne 0  and  Rmax.  The  strongest  dependence  is  on 
ne 0  [Denton  et  al ,  2002b].  Our  results  might  be  summarized 
by  two  simple  statements.  For  large  ne0,  such  as  occurs  in 
the  plasmasphere,  a  is  on  average  close  to  zero  (implying 
that  nc  is  roughly  constant  along  the  field  line).  When  ne0 
decreases,  so  does  nc  at  R  =  2  RE,  but  the  value  there  does 
not  decrease  below  ~100  cm"3.  This  is  apparently  because 
when  flux  tubes  become  depleted,  the  density  at  low  radius 
R  ~  2  RE  is  approximately  just  a  function  of  R,  independent 
of  Rmax.  (Note  that  the  field  lines  are  converging  at  low 
altitude  owing  to  the  geometry  of  the  dipole  field.)  Figure  9 
shows  a  simple  application  of  these  ideas.  The  thick  and 
thin  solid,  dashed,  and  dotted  curves  are  the  same  as  were 
plotted  in  Figure  7b.  The  dot-dashed  curves  in  Figure  9  are 
found  for  plasmasphere  (thick  dot-dashed  curve  in  Figure  9) 
and  plasmatrough  (thin  dot-dashed  curve  in  Figure  9)  from 


the  nc0  values  plotted  in  Figure  7a  in  the  following  way:  if 
ne 0  >  150  cm  ,  a  =  0,  while  if  «e0  <150  cm-3,  ne  at  R  =  2 
Re  is  150  cm-3,  and  a  is  found  from  equation  (3).  It  is  clear 
that  this  simple  model  (dot-dashed  curves  in  Figure  9) 
captures  the  general  trend  of  the  dependence  of  a  with 
respect  to  Rmax. 

5.  Kp  Dependence 

5.1.  Dependence  of  a 

[36]  As  discussed  in  section  1,  Denton  et  al  [2002b] 
found  that  they  could  model  a  using  equation  (2)  as  a 
function  of  «e0  and  i?max.  Here  we  examine  the  a  depen¬ 
dence  as  a  function  of  MLT  and  the  averaged  Kp,  (Kp). 
Following  Gallagher  et  al  [1988],  we  calculate  the  average 
at  the  current  time  t  by  averaging  over  earlier  times  /  using 
the  weighting  factor  exp(— (t  -  t')/t0),  where  t0  =  3.0  days. 
{Denton  et  al  [2002a]  used  1.5  days,  but  we  find  that  a 
3 -day  timescale  yields  a  parameter  that  correlates  better  to 
changes  in  density  as  described  in  sections  5.2  and  5.3.) 
Figure  10  indicates  that  there  is  no  clear  dependence  of  the 
average  a  —  amodei  on  MLT  or  (Kp).  (The  last  (Kp)  bin 
with  (Kp)  =  3.7  shows  a  -  amodei  slightly  higher  than  0, 
~0.5,  but  there  are  not  many  data  points  in  this  bin,  and  the 
scatter  in  a  -  amodei  at  this  value  of  (Kp)  is  very  large.) 
Since  there  is  a  significant  dependence  of  a  on  MLT  and 
(Kp)  (as  shown  by  Denton  et  al  [2002a]),  this  indicates  that 
the  MLT  and  Kp  a  dependence  can  be  accounted  for 
predominantly  by  changes  in  ne o- 

5.2.  Kp  Dependence  of  «e 

[37]  Figure  11  shows  the  (Kp)  dependence  of  the  electron 
density  nQ  at  the  magnetic  equator  (solid  curves)  and  at  R  = 
2  RE  (dashed  curves).  Figure  11a  shows  the  (Kp)  depen¬ 
dence  for  the  plasmasphere  data,  while  Figure  lib  shows 
the  (Kp)  dependence  for  the  plasmatrough  data.  The  three 


^max  «,) 


Figure  9.  The  thick  and  thin  solid,  dashed,  and  dotted 
curves  are  the  same  as  were  plotted  in  Figure  7b.  The  dot- 
dashed  curves  are  found  from  the  ne0  curves  in  Figure  7a 
(solid  curves  in  Figure  7a)  for  plasmasphere  (thick  dot- 
dashed  curve  in  this  figure)  and  plasmatrough  (thin  dot- 
dashed  curve  in  this  figure)  in  the  following  way:  if  «e0  > 
150  cm"3,  a  =  0,  while  if  ne 0  <150  cm"3,  «e  at  R  -  2  RE  is 
150  cm"3,  and  a  is  found  from  equation  (3). 
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Figure  10.  Plot  of  a  -  am(,rfcl  values  versus  (a)  MLT  and 
00  (Kp)  (average  of  Kp  with  a  timescale  of  3  days  as 
described  in  section  5.1)  for  all  categories  of  plasma  data 
combined.  The  solid  curves  show  the  average  values  in  bins 
(1  hour  bins  for  MLT),  while  the  dashed  line  in  Figure  10b) 
is  a  linear  fit. 


discussion  of  the  similarities  between  the  plasmatrough  and 
polar  cap.) 

[39]  The  smallest  spread  in  the  values  (width  between  the 
upper  and  lower  curves)  occurs  for  the  plasmatrough  data 
values  at  R  =  2  RE  (dashed  curves  in  Figure  lib).  This  is 
understandable  considering  that  the  plasmatrough  data  best 
represents  the  levelling  off  of  «e  at  large  Rmm  described  in 
section  4.2,  since  the  plasmatrough  data  more  typically  has 
a  low  value  of  ne  approaching  MOO  cm-3.  (On  the  other 
hand,  see  the  discussion  in  the  last  paragraph  of  section  4.2.) 

[40]  In  addition  to  the  (Kp)  dependence  shown  in 
Figure  11,  it  is  important  to  note  that  the  plasmasphere  data 
set  has  on  average  a  lower  value  of  (Kp)  =  1 .6  than  does  the 
plasmatrough  data  set,  with  (Kp)  =  2.3.  Noting  the  large 
difference  in  the  ne0  values  plotted  in  Figure  7a  for  the 
plasmasphere  (thick  solid  curve  in  Figure  7a)  and  plasma¬ 
trough  (thin  solid  curve  in  Figure  7a),  the  greatest  effect  of 


f 


curves  for  each  type  of  data  represent  the  log  average  value 
of  the  density  (middle  curve)  and  the  log  average  plus  or 
minus  one  standard  deviation  (upper  and  lower  curves). 
Thus  the  three  solid  curves  in  Figure  11a  show  the  log 
average  value  of  the  density  versus  (Kp)  for  the  plasma¬ 
sphere  at  R  =  Rmzx  (the  magnetic  equator)  and  the  log 
average  plus  or  minus  one  standard  deviation.  In  order  to 
reduce  the  effect  of  the  /?max  dependence,  we  have  fit  the  ne 
curves  in  Figure  7  with  functions  neL(Rmax),  and  the  density 
values  have  been  divided  by  this  function  (evaluated  at  the 
2^  max  value  corresponding  to  each  data  point)  before  binning 
the  nJneL  values  with  respect  to  (Kp).  Thus  the  overall 
vertical  position  of  each  set  of  three  nJneL  curves  in 
Figure  11  should  be  ignored.  The  important  thing  to  note 
is  the  (Kp)  dependence  of  each  set  of  three  curves  and  the 
relative  spread  of  the  three  curves. 

[38]  Figure  1  la  shows  that  for  the  plasmasphere,  the 
average  density  at  the  magnetic  equator  (middle  solid  curve) 
decreases  significantly  with  respect  to  (Kp)  (a  factor  of  2.0 
between  (Kp)  =  1 . 1  and  (Kp)  =  3.3).  This  is  to  be  expected  if 
low  (Kp)  indicates  that  more  refilling  has  been  allowed  to 
take  place.  The  average  density  at  R  =  2  RE  (middle  dashed 
curve)  shows  a  smaller  decrease  in  density  with  respect  to 
(Kp).  In  the  plasmatrough  the  average  density  at  the  mag¬ 
netic  equator  (middle  solid  curve  in  Figure  lib)  decreases  a 
small  amount  with  respect  to  (Kp)  (24%  between  (Kp)  =  1 .3 
and  (Kp)  =  3.5).  The  plasmatrough  density  at  R  =  2  RE 
(middle  dashed  curve  in  Figure  lib)  also  decreases  with 
respect  to  (Kp)  for  (Kp)  <  2.5  but  then  appears  to  increase 
with  respect  to  (Kp)  at  (Kp)  >  2.5.  Considering  the  large 
uncertainty  and  the  fact  that  there  are  only  six  events 
represented  in  the  two  data  bins  with  the  largest  values  of 
(Kp),  it  would  be  unwarranted  to  stress  this  result  too  much. 
However,  if  further  study  verifies  this  dependence,  the  trend 
at  low  (Kp)  could  be  caused  by  the  same  effect  seen  in  the 
plasmasphere  data,  that  at  lower  (Kp)  there  has  been  more 
refilling.  At  high  (Kp),  there  might  be  no  significant  refilling 
(due  to  the  large  activity),  and  in  that  case,  the  larger  density 
at  high  (Kp)  could  be  due  to  greater  ionospheric  outflow  at 
high  Kp.  A  study  of  the  polar  region,  where  the  magnetic 
field  lines  are  open,  shows  that  the  electron  density  is  higher 
at  large  Kp  [Nsumei  et  al.,  2003],  (See  section  3.1  for  a 


Independent  (or  rather  i?max-dependent)  function  as 
described  in  the  text.  The  solid  curves  represent  the  (Kp) 
dependence  for  the  density  at  R  =  /}max  (the  magnetic 
equator),  while  the  dashed  curves  give  the  (Kp)  dependence 
for  the  density  at  R  =  2  R%.  The  middle  curve  in  each  set  of 
three  (solid  or  dashed)  is  the  log  average  density  in  the  (Kp) 
bin,  while  the  upper  and  lower  curves  are  the  log  average 
plus  or  minus  one  standard  deviation. 
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(Kp)  on  the  density  is  likely  to  be  its  effect  on  the  type  of 
plasma  observed  (sphere  or  trough). 

5.3.  Kp  Dependence  of  «e  From  CRRES  Data 

[41]  While  Figure  11  shows  that  the  equatorial  density 
decreases  for  increasing  (Kp)  (particularly  for  the  plasma- 
sphere),  this  result  seems  to  be  contradicted  by  the  results  of 
Sheeley  et  al  [2001],  based  on  plasma  wave  measurements 
of  electron  density  observed  by  the  CRRES  spacecraft. 
Figure  6  of  their  paper  shows  that  the  electron  density  does 
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Figure  12.  Electron  density  ne  versus  Kp  for  the  plasma- 
sphere  using  the  CRRES  plasma  wave  data  set  [Sheeley  et  al , 
2001].  In  Figure  12a  the  log  average  electron  density  ne 
versus  (Kp)  (solid  curves)  and  versus  Kpmax  (dotted  curves) 
for  the  L  shell  ranges  indicated  in  the  figure.  In  Figure  12b  the 
log  average  of  values  of  ne  normalized  to  Sheeley  et  al.’s 
L-dependent  plasmasphere  model  density  (equation  (8)) 
versus  (Kp)  (dashed  curves).  The  middle  dashed  curve  is  the 
log  average  value  in  each  (Kp)  bin,  while  the  upper  and  lower 
dashed  curves  are  the  log  average  plus  or  minus  one  standard 
deviation.  The  solid  curves  are  the  same  as  the  dashed  curves, 
except  that  the  density  values  have  been  adjusted  to 
equatorial  values  ne0  using  equation  (1)  with  equation  (2). 


0  12  3  4  5  6 

Kp 


Figure  13.  Like  Figure  12,  except  for  the  plasmatrough 
data. 


not  depend  on  Kpmax,  where  Kpmax  is  the  maximum  value 
of  Kp  during  the  last  24  hours  before  each  density  mea¬ 
surement.  Because  of  this  difference  in  results,  we  have 
reanalyzed  the  CRRES  data  ourselves  (see  Sheeley  et  al 
[2001]  for  a  description  of  the  data).  Figures  12  and  13 
show  the  results  of  our  analysis  for  the  plasmasphere  and 
plasmatrough,  respectively.  Here,  we  follow  Sheeley  et  al. 
and  use  equation  (14)  as  the  borderline  density  separating 
plasmasphere  and  plasmatrough  data.  (As  discussed  in 
section  4.1,  this  borderline  value  works  fairly  well  for 
distinguishing  the  plasmasphere  and  plasmatrough  data; 
the  equation  (15  cm“3)  (6.6/Rmax)4  does  even  better.)  The 
dotted  curves  in  Figure  12a  show  the  log  average  electron 
density  ne  versus  Kpmax  for  the  L-shell  ranges  indicated  in 
the  figure.  These  curves  should  be  similar  to  the  curves  in 
Figure  6b  of  Sheeley  et  al  [2001].  However,  unlike  Sheeley 
et  al.,  we  find  a  clear  trend  of  decreasing  ne  with  respect  to 
Kpmax,  especially  at  the  lower  L  values.  The  solid  curves  in 
Figure  12a  show  the  log  average  nc  versus  (Kp)  (averaged 
with  a  3-day  timescale;  see  section  5.1).  It  is  clear  that  there 
is  even  a  stronger  dependence  of  nQ  on  (Kp)  than  on  Kpmax, 
and  therefore  (Kp)  is  a  better  indicator  of  the  Kp  depen¬ 
dence.  This  is  not  suiprising  considering  that  (Kp)  is  a 
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measure  of  the  cumulative  activity  as  opposed  to  the  activity 
during  one  particular  hour  of  time. 

[42]  The  solid  curves  in  Figure  1 2b  are  roughly  equivalent 
to  the  solid  curves  in  Figure  11.  Here,  the  ne  values  have 
been  normalized  to  Sheeley  et  al.’s  Z-dependent  plasma- 
sphere  model  density  (equation  (8)),  and  the  values  of  nJneL 
are  plotted  versus  (Kp).  Since  CRRES  did  not  have  a  polar 
orbit,  we  do  not  have  paired  measurements  at  one  value  of  L 
with  different  radii  (CRRES  had  an  orbit  inclined  18°  from 
file  geographic  equator).  Therefore  in  order  to  calculate  the 
equatorial  density,  we  have  used  equation  (1)  with  equation 
(2).  (In  this  subsection,  the  dipole  field  model  is  used.)  The 
middle  solid  curve  is  the  log  average  value  of  nJneL,  while 
the  upper  and  lower  curves  are  the  log  average  value  plus  or 
minus  one  standard  deviation.  The  dashed  curves  are  the 
same  except  that  the  measured  «e  values  are  used  with  no 
correction  to  ne q.  In  this  case,  the  solid  and  dashed  curves 
are  almost  indistinguishable.  This  is  because  at  the  high 
densities  characteristic  of  the  plasmasphere,  a  is  small  so 
that  equation  (1)  leads  to  roughly  constant  density  along  the 
field  line  and  the  difference  between  the  locally  measured  nc 
and  ne 0  is  small.  (The  unadjusted  ne  values  were  also  used  in 
Figure  12a.)  The  decrease  in  nJneL  with  respect  to  (Kp)  in 
Figure  12b  is  larger  than  the  uncertainty,  which  helps  us 
have  some  confidence  that  the  Kp  dependence  for  the 
plasmasphere  density  is  real.  Over  the  same  range  of  (Kp) 
as  was  plotted  in  Figure  II  (1.1— 3.3),  nJneL  decreases  by  a 
factor  of  1.6,  somewhat  smaller  than  the  decrease  in 
Figure  1 1  (a  factor  of  2.0). 

[43]  Figure  13  shows  the  same  quantities  as  were  plotted 
in  Figure  12,  except  for  the  CRRES  plasmatrough  data  set 
(ne  values  lower  than  the  value  of  nc  given  by  equation 
(14)).  From  the  dotted  curves  in  Figure  13a,  we  can  see  that 
there  is  no  apparent  dependence  of  ne  versus  Kpmox.  Neither 
is  the  dependence  of  ne  on  (Kp)  very  clear  (solid  curves  in 
Figure  13a).  On  the  other  hand,  when  the  density  values  are 
normalized  to  Sheeley  et  al.’s  [2001]  plasmatrough  model 
including  MLT  dependence  (their  equation  (7)), 

»eo  =  124  cm-J(3/Rmax)40, 

+  36  cm  3(3 /Rmn)  cos  i^MLT  —  j7.7(3//tmax)2'°+12j  j 

W12}),  (16) 

there  is  an  apparent  dependence  of  the  normalized  nj 
»er-MLTon  (Kp)  (Figure  13b;  here  nei_MLT  is  dependent  on 
both  L  and  MLT).  The  biggest  factor  causing  nJneL_M LT  to 
decrease  with  respect  to  (Kp)  is  the  MLT  dependence;  if  the 
MLT-dependent  cosine  term  in  equation  (16)  is  neglected, 
nJneL- mlt  does  not  decrease  with  respect  to  (Kp).  The 
MLT  dependence  makes  a  difference  in  the  results  because 
the  average  Kp  values  were  higher  during  the  middle  of 
1991  when  CRRES  was  sampling  the  afternoon  MLT  values 
than  when  it  sampled  dawn  and  nighttime  MLT  values.  In 
fact,  the  average  value  of  Kp  was  2.2  during  the  first  part  of 
the  CRRES  mission,  1  August  1990  to  23  March  1991,  but 
3.2  from  the  24  March  1991  storm  to  the  end  of  the  mission 
on  12  October  1991.  CRRES  did  not  completely  sample  all 
MLT  values  during  its  entire  mission;  the  apogee  of  the  orbit 
started  in  the  dawn  MLT  sector  and  moved  through  midnight 
to  afternoon  MLT  values  at  the  end  of  the  mission.  The  large 
density  values  at  (Kp)  ~  4.5  in  Figure  13a  are  predominantly 


at  afternoon  values  of  MLT,  where  K£z,_mlt  (equation  (16)) 
has  a  peak  with  respect  to  MLT.  (This  correlation  between 
Kp  and  MLT  for  the  CRRES  data  will  obviously  complicate 
modeling  both  the  Kp  and  MLT  dependence.) 

[44]  As  seen  from  Figure  13b,  the  total  decrease  in  nj 
«e£— MLT  over  the  range  of  (Kp)  plotted  is  a  factor  of  1.6  (as 
compared  to  2.2  for  the  plasmasphere).  Note  also  that  for 
the  plasmatrough  data,  the  use  of  equation  (1)  with  equation 
(2)  to  find  «e0  makes  a  greater  difference  than  it  did  for  the 
plasmasphere,  as  indicated  by  the  difference  between  the 
solid  and  dashed  curves  in  Figure  13b.  This  is  because 
the  values  of  neo  are  lower  for  the  trough,  so  ot  from 
equation  (2)  is  larger. 

[45]  The  plasmatrough  density  may  be  dependent  on  Kp 
as  indicated  by  Figure  lib  and  Figure  13b,  but  considering 
the  small  amount  of  data  at  large  Kp  represented  in 
Figure  lib  and  the  lack  of  any  Kp  dependence  in 
Figure  13a,  we  are  hesitant  to  come  to  a  firm  conclusion 
about  this  dependence  without  further  study.  However,  the 
dependence  in  the  plasmasphere  is  more  clear.  We  would 
most  expect  to  see  Kp  dependence  for  the  plasmasphere, 
since  the  plasmasphere  develops  from  a  period  of  refilling. 
Furthermore,  the  observed  Kp  dependence  can  account  for 
the  difference  between  the  average  plasmasphere  model 
from  our  results  (thick  solid  curve  in  Figure  6,  which  we 
have  fit  by  equation  (5)  with  equation  (9))  and  that  of 
Sheeley  et  al.  [2001]  (thick  dashed  curve  in  Figure  6,  found 
from  equation  (8)).  The  average  value  of  (Kp)  for  CRRES’s 
plasmasphere  data  (density  measurements  with  ne  greater 
than  the  value  of  ne  given  by  equation  (14))  is  2.5,  whereas 
the  average  value  for  our  Polar  plasmasphere  data  set  is  1.6. 
On  the  basis  of  Figure  11,  this  would  cause  the  CRRES 
average  density  to  be  -30%  lower  than  the  Polar  average 
density,  and  this  is  about  the  difference  between  the  two 
curves  in  Figure  6  (thick  solid  and  thick  dashed  curves).  If 
the  average  value  of  (Kp)  is  even  lower  for  the  plasma¬ 
sphere  data  of  Carpenter  and  Anderson  [1992],  the  fact  that 
their  density  is  higher  than  the  other  models  (at  least  for 
R-mm  >  4  Re;  see  Figure  6)  could  be  accounted  for. 

6.  Summary 

[46]  The  Polar  PWI  data  set  spans  the  parameters  2.5  Re  < 
Rmm  ^  8.5  Re,2Re<  R  <  Rmax,  and  2  <ne0<  1500  cm-3, 
and  this  data  set  has  been  used  to  generate  a  model  for  the 
average  c*  (equation  (2))  as  a  function  of  equatorial  density 
ffeo  and  Rm ax  assuming  the  power  law  form  for  ne  (equation 
(1))  [Denton  et  al.,  2002b],  Variations  of  the  average  a  with 
respect  to  MLT  or  Kp  [, Denton  et  al.,  2002a]  can  be  attributed 
predominantly  to  changes  in  n^  (Figure  10).  Typical  values 
of  a  were  0-1  in  the  plasmasphere  and  2-3  in  the  plasma- 
trough  (Figure  7b).  While  the  low  a  values  within  the 
plasmasphere  are  consistent  with  diffusive  equilibrium,  the 
larger  a  values  characteristic  of  the  plasmatrough  have  a 
dependence  which  is  intermediate  between  diffusive  equi¬ 
librium  and  a  collisionless  model  (roughly  ne  oc  R~4) 
[Lemaire  and  Gringauz,  1998],  It  is  important  to  note  that 
oui  results  describe  only  the  average  field  line  dependence; 
there  will  surely  be  deviations  from  the  average  behavior! 
some  of  which  will  be  related  to  time-dependent  effects. 

[47]  We  have  examined  several  events  with  electron 
density  «e  measured  on  large  and  small  radius  portions  of 
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an  orbit  (Figures  2-5).  These  events  exhibit  features  that  we 
commonly  find  in  the  Polar  data.  In  particular,  the  plasma- 
pause  is  much  more  evident  on  the  part  of  the  orbit  with 
large  radius  (near  the  magnetic  equator)  than  on  the  part  of 
the  orbit  with  small  radius  (R  =  2-3  RG),  and  on  the  small 
radius  portion  of  the  orbit  where  R  ~  2  RG,  ne  often  levels 
off  to  nG  ~  100  cm-3  at  large  Rmax.  This  is  consistent  with 
the  results  of  Reinisch  et  ah  [2004]  for  plasmaspheric  field 
lines,  who  show  that  temporal  variations  associated  with 
flux  tube  depletion  and  refilling  are  less  at  low  altitude  than 
at  the  magnetic  equator.  We  also  present  an  example  of 
plasmaspheric  plasma  extending  out  to  at  least  L  ~  9  on  the 
dawnside  during  particularly  calm  geomagnetic  conditions 
(as  indicated  by  Kp;  see  section  3.3). 

[48]  We  have  calculated  the  average  equatorial  profile  of 
«e o  versus  Rmax  for  the  plasmasphere  data  (including  the 
“gradually  decreasing”  category  of  Denton  et  ah  [2002a] 
(Appendix  A)  and  find  that  it  lies  between  the  plasmasphere 
models  of  Carpenter  and  Anderson  [1992]  and  Sheeley  et  al. 
[2001].  Our  average  plasmatrough  density  scales  with  re¬ 
spect  to  Rmax  like  R~ 3X  .  However,  excluding  the  largest  Rmax 
values  >6  RG  (at  which  we  have  reason  to  doubt  the  accuracy 
of  our  average  values;  see  section  4.1),  our  plasmatrough  «e0 
agrees  well  with  Sheeley  et  al ’s  [2001]  result  R~ ax°  or  with 
Carpenter  and  Anderson’s  R^ix  (Figure  6). 

[49]  In  addition  to  the  average  equatorial  profile,  we  also 

calculated  the  log  average  nG  values  atR~2Re  versus  Rmax 
(dashed  curves  in  Figure  7a).  Like  the  density  profiles  for 
the  example  events,  these  also  level  off  at  ~  100  cm’3 

at  large  Rmax.  We  also  found  profiles  of  nG  at  R  =  2  Re  by 
using  equation  (2)  with  our  average  nG0  profiles  as  an  input 
(dotted  curves  in  Figure  7a).  These  also  approached  values 
of  nG  ~  100-200  cm’3  at  large  Rmax.  The  same  result  was 
also  observed  in  the  meridional  structure  of  our  plasma¬ 
trough  results  (equation  (12)  with  equation  (2)),  as  shown  in 
Figure  8.  These  properties  are  similar  to  those  of  the  density 
in  the  polar  cap,  where  the  density  is  mostly  a  function  of 
radius  [Nsumei  et  al. ,  2003].  However,  the  field  line 
dependence  and  Kp  dependence  for  the  plasmatrough  are 
not  the  same  as  in  the  polar  cap  (sections  3.3  and  5.3). 

[so]  There  were  some  differences  in  the  ne  values  at  R  ~ 

2  RG,  depending  on  how  they  were  calculated  (log  average 
values  or  the  value  found  using  equation  (2)  with  «e0  as  an 
input).  Correspondingly,  there  were  also  differences  in  the  a 
values  computed  in  different  ways  (average  value,  value 
from  equation  (2),  or  value  found  using  equation  (3)  with 
the  average  nG 0  and  value  of  nG  for  R  =  2  RG;  these  are 
plotted  in  Figure  7b).  These  differences  demonstrate  the 
limitations  of  our  power  law  model. 

[51]  We  need  to  keep  in  mind  the  fact  that  some  low- 
density  measurements  at  low  radius  could  be  missed  be¬ 
cause  fp!fce  <  1  (section  3.4).  This  could  bias  our  results 
somewhat.  In  particular,  the  density  at  low  radius  and  large 
i?max  might  have  in  some  cases  decreased  significantly 
below  100  cm-3,  and  PWI  could  not  measure  it.  On 
the  other  hand,  we  often  see  the  density  level  off  at  nG 
~  100  cm-3  at  small  radius  and  large  Rmax  when  the 
equatorial  density  is  still  decreasing  with  respect  to  Rmax 
(Figures  2a  and  5a),  and  nG  ~  100  cm’3  at  R  2  Re  is 
consistent  with  the  average  field  line  dependence  using 
density  measurements  at  higher  altitude  (Figure  7  and 
discussion  in  the  last  paragraph  of  section  4.2). 


[52]  We  also  presented  a  simplified  model  for  a  using  the 
observations  that  at  large  equatorial  density  ne0,  the  density 
is  roughly  constant  along  field  lines  and  that  nG  at  R  ~  2  RG 
approaches  as  a  rough  lower  limit  100  cm’3.  This  simple 
model  (dot-dashed  curves  in  Figure  9)  captures  the  general 
trend  of  the  dependence  of  a  with  respect  to  RmaX}  though  the 
agreement  with  the  other  values  of  a  is  certainly  not  exact. 

[53]  Finally,  we  examined  the  Kp  dependence  of  the 
equatorial  density  and  density  at  R  ~  2  Re.  For  the 
equatorial  density,  we  supplemented  our  Polar  data  set  with 
data  from  the  CRRES  spacecraft  [Sheeley  et  al,  2001].  The 
plasmasphere  nG  values  at  the  magnetic  equator  clearly 
decrease  with  respect  to  (Kp)  (Kp  averaged  with  a  3-day 
timescale  as  described  in  section  5. 1),  as  shown  in  Figures  1 1  a 
and  12.  We  also  showed  that  the  quantity  (Kp)  better 
indicates  the  variation  of  density  with  respect  to  Kp  than 
the  maximum  value  of  Kp  during  the  last  24  hours 
(section  5.3).  The  fact  that  the  plasmasphere  density 
decreases  with  respect  to  (Kp)  is  consistent  with  the  idea 
that  greater  refilling  occurs  over  a  period  of  low  geomag¬ 
netic  activity.  This  dependence  of  the  plasmasphere  density 
on  (Kp)  can  also  explain  the  difference  in  the  average 
density  based  on  the  Polar  data  and  that  of  Sheeley  et  al.,  as 
shown  in  the  last  paragraph  of  section  5.3.  The  dependence 
of  the  plasmatrough  density  with  respect  to  (Kp)  is  not  so 
clear  and  perhaps  more  complicated. 

Appendix  A:  Definitions  of  Plasmasphere  and 
Plasmatrough  Data 

[54]  Denton  et  al  [2002b]  separated  the  Polar  PWI  data 
into  three  categories.  If  there  was  a  plasmapause  (drop  in  ne  of 
at  least  a  factor  of  3  within  AKmax  =  0.4),  they  categorized  the 
density  data  inside  (outside)  the  plasmapause  as  plasma¬ 
sphere  (plasmatrough)  data.  If  the  orbit  had  only  ne  values 
>300  cm-3,  they  categorized  it  as  plasmasphere  data.  The 
plasmatrough  data  was  almost  exclusively  identified  as  the 
region  outside  a  plasmapause;  however,  in  some  cases  no 
plasmapause  was  discemable,  but  the  density  was  so  low  that 
the  data  was  categorized  as  plasmatrough  data.  For  instance, 
in  one  case,  ne  was  M0  cm-3  at  Rmax  ~4RG  and  decreased 
gradually  down  to  values  ~  1 0  cm’3  at  larger  Rmax.  Sometimes 
ne  decreased  gradually  with  respect  to  RmSLX  from  values 
>300  cm’3  to  values  ~10  cm’3  with  no  clear  plasmapause 
(Figure  5),  and  they  counted  this  data  as  a  third  “gradually 
decreasing”  (no  plasmapause)  category.  In  most  cases  in  the 
gradually  decreasing  category,  the  density  at  low  Rmax  was 
higher  than  300  cm’3,  M000  cm’3  at  Rmax  ~  3  RG  and 
MOO  cm”3  at  Rmax  ~  3  RG.  For  the  purposes  of  this  paper, 
the  plasmasphere  and  “gradually  decreasing”  categories  are 
combined  into  a  new  plasmasphere  category.  The  “gradually 
decreasing”  plasma  is  assumed  to  be  extended  plasmasphere. 

Appendix  B:  Details  Concerning  Figure  5 

[55]  Profiles  of  density  with  respect  to  Rmax  similar  to  that 
shown  in  Figure  5  were  observed  at  the  same  MLT  ~  7- 
8  hours  on  the  next  two  Polar  passes,  on  20  December  1996 
at  1615-1826  UTand  on  21  December  at  0950-1210  UT. 
Similar  profiles  of  density  versus  Rmax  were  also  observed  at 
dusk  MLT  ~  20  hours,  before  and  after  the  time 
corresponding  to  the  event  in  Figure  5  (at  19  December, 
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0747-1010  UT,  and  at  20  December,  0100-0310  UT,  20 
December,  1853—2045  UT,  and  21  December,  1248- 
1435  UT).  The  20  December,  0100-0310  UT  observation 
included  “gradually  decreasing”  density  out  to  at  least 
•^max  ~  12  Re.  On  the  other  hand,  a  plasmapause  at  i?max  = 
5  Re  with  a  drop  in  density  of  a  factor  of  5  was  observed  at 
MLT  ~  0800  (same  MLT  as  Figure  5)  at  17  hours  (UT) 
earlier  than  the  time  corresponding  to  Figure  5.  In  order  to 
inteipret  the  “gradually  decreasing”  condition  (no  plasma- 
pause)  as  applying  uniformly  to  all  azimuths,  it  is  necessary 
to  assume  that  the  density  increased  significantly  in  the 
region  just  outside  Rmax  =  5  Re  from  ~20  cm-3  to 
~200  cm  3  during  a  period  of  17  hours.  The  existence  of 
this  extended  plasmasphere  density  correlated  with  an 
extended  period  of  low  geomagnetic  activity  as  indicated 
by  a  low  value  of  Kp.  Data  from  NSSDC  OMNIWeb 
(NOAA)  shows  that  the  hourly  average  Kp  on  19  December 
1996  was  equal  to  0.7  from  0000  to  0200  UT,  0  from  0300  to 
2000  UT,  0.7  from  2100  to  2300  UT,  and  0  again  at  2400  UT. 
At  2300  UT,  the  value  of  Kp  averaged  using  the  technique  of 
section  5.1  (with  a  timescale  of  3.0  days),  (Kp),  was  1.0. 

Appendix  C:  Total  Flux  lYibe  Particle  Content 

[56]  The  total  flux  tube  particle  content  (per  flux)  will  be 
constant  across  the  plasmatrough  under  the  assumptions  that 
(1)  the  plasmatrough  flux  tubes  are  simultaneously  depleted 
at  some  time,  (2)  they  subsequently  begin  to  refill  due  to 
upward  flux  of  particles  from  the  ionosphere,  (3)  a  similar 
fraction  of  that  upward  flux  is  trapped  in  each  flux  tube,  and 
(4)  the  plasmatrough  flux  tubes  map  down  to  a  similar 
latitude  at  the  Earth’s  surface.  The  last  condition  results 
from  the  geometry  of  the  nearly  dipolar  magnetic  field. 
Because  the  flux  tubes  map  down  to  a  similar  latitude  at  the 
Earth’s  surface,  B  at  the  Earth’s  surface  is  roughly  constant 
with  respect  to  L  so  that  the  area  of  a  flux  tube  (per  magnetic 
flux)  at  the  Earth’s  surface  is  roughly  constant.  The  other 
consequence  of  the  flux  tubes  mapping  to  a  similar  latitude 
is  that  the  upward  refilling  flux  of  particles  (particles  per 
area  per  time)  from  the  ionosphere  will  be  roughly  constant 
with  respect  to  L.  Then  the  particles  per  time  filling  the  flux 
tube  will  be  constant  with  respect  to  L. 
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